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Abstract. Pulse width modulation (PWM) inverter drive system has been widely 
applied for its great speed-control characteristics. However, through coupling 
capacitances in AC motors, the high-frequency common-mode voltage caused by 
the inverter will induce bearing currents, which leads to bearing electro erosion and 
eventually premature failure. Therefore, it is necessary to investigate the mech- 
anism of bearing currents. On this basis, the influence of the bearing currents is 
predicted and corresponding mitigation techniques are proposed. This article intro- 
duces the modeling process of typical bearing currents, including capacitive dv/dt, 
EDM, circulating, and rotor grounding types. The calculations of coupling capac- 
itances in AC motors and the bearing capacitance are also presented. Furthermore, 
the commonly used mitigation techniques for bearing currents are summarized. 
Finally, the future research direction of the bearing currents is prospected. 
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1 Introduction 


In line-operated electrical machines, unbalanced magnetic fields caused by rotor eccen- 
tricities and axial cooling holes in stator and rotor laminations can create an alternating 
magnetic flux encircling the rotor shaft, which further induces alternating voltages. Since 
the loop impedance is relatively small, a circulating current is generated in the loop “‘sta- 
tor frame-drive-end bearing-shaft-non-drive-end bearing”. Such currents are referred 
to as classical bearing currents [1-3]. With the increase of motor size, these classical 
bearing currents are more likely to occur due to the increased alternating magnetic flux 
and decreased lubrication impedance, which may lead to premature bearing failures. In 
the past decades, systematic researches on classical shaft currents have been carried out 
by many scholars and bearing companies, and standard limits have been given for the 
shaft voltage at line-operation [4]. With the improvement of motor manufacturing details 
and the application of mitigation techniques such as insulated bearings and grounding 
carbon brushes, the classical bearing current problem has been basically solved. 
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Due to the rapid development of high-performance semiconductor switching devices 
and advanced control strategies, AC motors operated with pulse width modulation 
(PWM) are widely applied to variable speed drive systems for their reduced torque rip- 
ples, increased system dynamics, and high efficiency. However, fast switching inverter 
operated motors are submitted to additional bearing currents that had not been of influ- 
ence before [5, 6], especially in the field of metro vehicle traction systems [7—9]. In 
recent years, statistics show that 40% of motor failures are caused by motor bearing dam- 
age, and 25% of motor bearing damage is caused by inverter-induced bearing currents 
[10]. Generated by the PWM operation of the inverter, the high-frequency common- 
mode (CM) voltage acts on the capacitive coupling path provided by three main stray 
capacitances in the motor (the stator winding-to-frame capacitance Cyr, rotor-to-frame 
capacitance C,f, and stator winding-to-rotor capacitance Cwr), causing the bearing volt- 
age vp. When vp surpasses the threshold voltage of the oil film, an electrical discharge 
will occur and the local high temperature may melt the metal causing damages such 
as microscopic pits and grooves, which is the main cause of bearing degradation and 
failure [11]. Thus, the phenomena of the inverter-induced bearing currents should be 
comprehensively investigated. 

So far, the researches on bearing currents in variable speed drive systems mainly 
focus on four aspects, including their generation mechanism, equivalent circuit model, 
related capacitances computation, and mitigation techniques. In this paper, an overview 
of the bearing currents phenomena has been presented, recent researches of the above 
aspects have been summarized, and the future research direction of the bearing currents 
has been discussed. 


2 Generation Mechanism of Bearing Currents 


S. Chen, D. Busse, and J. Erdman [12-14] first paid attention to the phenomena of 
bearing currents in variable speed drive systems. They summarized that the non-zero 
high-frequency CM voltage generated by the PWM inverter is the source of bearing 
currents. Chen [12] verified the major parasitic coupling path existing in the motor, 
which couples through parasitic capacitances from stator windings to the rotor body 
and then returns through the motor bearings to the commonly grounded stator frame 
as a closed-loop circuit. On this basis, Busse and Erdman [13, 14] pointed out that the 
bearing acts as a capacitor at rated motor speed due to hydrodynamic effects. In this 
case, the bearing voltage mirrors the CM voltage by a capacitive voltage divider and 
small capacitive dv/dt bearing currents appear at each switching instant of the inverter. 
However, if the threshold voltage of the lubricating film is surpassed, the film discharges, 
causing harmful electrostatic discharge (EDM) currents. 

Later, Annette Muetze and Andreas Binder published a series of papers [15-18] that 
systematically summarize the generation mechanism of dv/dt and EDM bearing currents 
in inverter-supply motors. As is stated in references, the capacitive dv/dt currents (5— 
200 mA) are considered to be harmless since they are so small, whereas small inverter-fed 
AC motors of up to typically 20 kW are likely to suffer from EDM bearing currents (0.5- 
3 A). Further, in [17, 19, 20], Muetze focused on the high-frequency circulating bearing 
currents in larger motors. The parasitic capacitance Cwis excited by the high dv/dt at the 
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motor terminals, creating a high-frequency ground current, which generates a circular 
flux around the motor shaft, inducing the shaft voltage. If the voltage is large enough 
to destroy the insulating properties of the bearing lubricating film, a circulating bearing 
current (0.5—20 A) along the loop of “stator frame- non drive end- shaft- drive end” 
appears, which mirrors the ground current and is of opposite direction in both bearings. 

In some configurations, the rotor is connected to the ground with a lower impedance 
path than the grounding of the stator frame [15, 16]. Thus, the bearing voltage is recon- 
figured, which results in a considerable rotor ground current flowing through the motor 
bearings. These currents can reach 1-35 A and prematurely damage the bearings. 

Based on the above references, the bearing currents can be divided into four categories 
according to their generation mechanisms, and their coupling paths are shown in Fig. 1 
[21]. 


frame 


stator iron 


rotor ground 


bearing currents 


Fig. 1. Schematic diagram of coupling paths of bearing currents. 


3 Equivalent Circuit Model of Bearing Currents 


3.1 dv/dt and EDM Bearing Currents 


The equivalent circuit models of dv/dt and EDM bearing currents are first proposed 
by Chen [22] and Busse [14]. Based on the transmission line theory, the three-phase 
distributed parameter circuit is eventually equivalent to the single-phase lumped param- 
eter circuit in [22]. However, the physical meanings of the parameters in this model are 
ambiguous, so it needs to be fitted according to the measured bearing current wave- 
forms. In [14], a zero sequence equivalent circuit model is proposed, as is shown in 
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Fig. 2. Moreover, the relationship between CM voltage and stator neutral to the ground 
is approximated by a second-order RLC circuit model, which can be expressed as (1): 


Vb Cwr 


BVR = = : 
Vng Cwr + Crp + 2Cp 


(1) 


Z series 


Fig. 2. Equivalent circuit model for dv/dt and EDM bearing currents. 


Figure 2 allows for the investigation of CM chokes or transformers, line reactors, 
and long cables through the modification of the series and parallel impedance elements 
(Zseries and Zparalle!), and the zero-sequence impedance of the stator windings is repre- 
sented by Zo. Although greatly simplified, the equivalent circuit provides a useful tool 
for the analysis of system parameters and their effect on bearing voltages without the 
complexity of the distributed system. 

On this basis, Annette Muetze [15, 23] proposed a simplified model of EDM currents 
and used a controllable switch K in Fig. 2 to simulate the oil film breakdown process. 
According to [24, 25], the maximum amplitude of EDM currents can be evaluated by 
the breakdown voltage and Rp. Based on this model, Muetze analyzed the error between 
the measured and calculated bearing currents and discussed the influence of different 
system parameters. In reference [26, 27], Oliver Magdun combined the high-frequency 
model of induction motor with the above-mentioned bearing current model and obtained 
the circuit model parameters conveniently and accurately by measuring the CM and 
differential mode impedance curves of the system. 


3.2 Circulating Bearing Currents 


Based on the generating mechanism analysis in [28], the equivalent circuit models of 
circulating bearing currents are first proposed by Maki-ontto and Luomi [29, 30]. With 
the high-frequency common-mode current icm causing the common-mode flux ® eir that 
varies with time, thereby inducing a voltage vsh, and with the circulating bearing current 
ip_cir flowing as a result of this voltage, the machine can be considered as a current 
transformer. Thus the bearing current ratio can be expressed as [19, 20]: 


lb_ cir 


BCR = 


(2) 
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where Lg is the equivalent magnetizing inductance, Lp_Fe, Rp_Fe, Lp_air are the equivalent 
impedances of the bearing current path. 


3.3 Rotor Ground Bearing Currents 


The equivalent circuit models of rotor ground bearing currents are first proposed by 
Guttowski [31]. In order to study the influence brought by the existence of rotor ground- 
ing path, the stator grounding impedance Zfg and rotor grounding impedance Zrg are 
added on the basis of the aforementioned EDM bearing current model. In [32], the effect 
of stator and rotor grounding impedances on the bearing voltage is studied through a 
Wheatstone bridge theory-based equivalent circuit. However, this model cannot consider 
the influence of rotor grounding impedance on the stator neutral to ground voltage. So 
far, most studies on rotor ground bearing currents remain in the qualitative analysis due 
to the complexity of the high order equivalent model. 


4 Related Capacitances Computation of Bearing Currents 


If the parameters in the equivalent circuit can be theoretically calculated according to 
the motor and bearing structures, it will bring great convenience to the prediction of the 
bearing voltages and currents. Numerous researches have focused on the theoretical cal- 
culation of parasitic capacitance and bearing capacitance parameters in non-circulating 
bearing current models. In contrast, the circulating bearing current model and its param- 
eter calculation method are relatively unified, which are systematically introduced in 
literature [19] and will not be described in this paper. 


4.1 Motor Capacitances 


Theoretical calculation methods of motor capacitances mainly include the analytical 
method [14, 15, 33, 34] and the finite element method [34-36]. Figure 3 shows the 
simplified schematic diagram of motor high-frequency equivalent capacitance used in 
the analytical method. 


slot 


frame 
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stator tooth 
air gap (om 
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Fig. 3. Simplified schematic diagram of motor high-frequency equivalent capacitances 
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The analytical calculation of Cwr is first introduced in [14]. Cwf is modeled as a 
number of Ns parallel plate capacitors. Each plate capacitor is regarded as a rectangular 
conductor of length lfe, width Wq and depth W,, the dielectric between stator winding 
and frame is €;;, and the distance is dsjot. In [15], a form factor Fe is introduced to take 
into account the irregular surface of the windings. Thus, the analytical expression of Cwf 
is derived as 


Cwt = FeNséri€o(Wa + 2Ws)lFe/dstot (3) 


where €,9 = 8.854-e*(— 12) [F/m] is the permittivity of vacuum. 

In addition, the parasitic capacitance of stator windings-to-stator teeth is considered 
in [37], which makes the results more accurate. Different from the above literature, 
Magdun [26] started from the perspective of high-frequency common-mode impedance 
characteristics of induction motors. Below the first anti-resonance frequency, the motor 
impedance characteristics were mainly determined by Cyr. This method does not involve 
simplification of the model and is thus more accurate. 

The analytical calculation of Cy is first introduced in [14]. The capacitance Cy 
is assumed to form a coaxial air-insulated capacitor between the stator iron and rotor 
surfaces with the length of stator stack lre, 


Crt = 270 €olFe/In (dsi/ dre) (4) 


where dsi is the stator inner diameter and dre is the rotor outer diameter. 

In [15], the effect of the stator slot opening is considered by Carter coefficient ke. 
Since the air gap is relatively small compared to the rotor diameter, the capacitor can be 
simplified as an air-insulated plate capacitor. 


Crt = £olfen dre / (koô) (5) 


The capacitance Cwr is formed between the stator windings in the slots and the rotor 
surface, and its analytical calculation method is first introduced in [14]. It assumes that 
Ns stator slots form parallel plate capacitors with the rotor surface, with a plate width 
of W, and length lfe. The distance between the plates is ô. 


Cwr = Ns€0WrlFe/6 (6) 


Obviously, the formula fails to consider the influence of stator slot opening and slot 
wedge. Subsequently, some improvements were made in the literature [15]. The distance 
between the plates comprises the stator tooth tip thickness ho, the air gap ô, the distance 
from the stator-slot to the winding dsjot, and the slot insulation thickness hy. Thus, a 
series connection of capacitors is formed as follows. 


Ns€0bolfe _ Ns&r1€0bolfe 1 


C. = —s rq = 5 eS — 
me ii hw + dslot a 1/Cwro + 1/Cwrt 


ô + ho 


where bg is the stator tooth width. 
On this basis, the capacitance of stator windings-to-stator teeth-to-rotor surface is 
taken into account in [37]. 


678 Y. Zhi et al. 


Magdun [36] first employed the 2D finite element method to calculate motor high- 
frequency capacitances and found that the calculation of Cyr and Cwr by the solid 
conductor model is not accurate enough. Therefore, scattered line model and 3D finite 
element method are proposed to improve the calculation accuracy [34-36]. 


4.2 Bearing Capacitance 


The theoretical calculation of bearing capacitance is first proposed in [14]. The model 
assumes Np pairs of spheres, and each pair includes an inner ball sphere within an outer 
raceway sphere. The radius of the roller element is Rpay, the distance is bg, and the 
relative permittivity of lubrication is €,2. 


Cy = No4 £0812/(1/Rpan — 1/Roa + bg) (8) 


In the above method, the improper equivalent of the sphere capacitor will lead to a 
large error of capacitance value. Subsequently, Muetze [15] used the minimum thickness 
of the lubrication film of the bearing Ayertz and the Hertz’ian contact area Syertz that is 
given by the elastic deformation of the balls of a roller bearing under mechanical pressure 
in practical operating conditions to calculate the bearing capacitance (Fig. 4). 


Fig. 4. Hertz contact model of bearing capacitance 


Cb = 0.5¢,€0SHertz / Hertz (9) 
Y > outer race Se | 
OY LLB : | : 
Cur Suertsz | 
sa | hertz ! 
Cre > >i 


inner race 


However, the corresponding electrode area calculated by the Hertz’ian contact area 
Shertz 18 smaller than the actual area, which leads to a smaller Cy. As a supplementary, 
Magdun [25] and Liu [38] considered the electrical coupling air capacitance Cair between 
the rolling elements and the raceway, which makes the calculation much more accurate. 

In addition, a new calculation method on the bearing capacitance is put forward in 
[39]. Based on the central film thickness of bearing oil obtained through the elastohy- 
drodynamic lubrication analysis, the bearing capacitance can be accurately calculated 
with the help of ANSYS Maxwell finite element analysis software. 


5 Mitigation Techniques for Bearing Currents 


According to the generation mechanism and equivalent circuit model of bearing currents, 
different mitigation techniques have been proposed to eliminate bearing currents. These 
techniques can be divided into two groups depending on their function. 
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The first group consists of mitigation techniques on the inverter side, which aim at 
the inverter-induced high-frequency CM voltage. The most straightforward method is to 
reduce the switching frequency at the cost of deteriorated control performance [40]. In 
[41, 42], the CM voltage is successfully suppressed by novel PWM control strategies, 
which are realized by software without additional hardware costs. From the perspec- 
tive of inverter topologies, multilevel inverters can naturally achieve the suppression of 
CM voltage compared with traditional two-level inverters [43, 44]. In addition, output 
filters intended for bearing voltage elimination are designed in [45-47]. Considering 
the complex design procedure of most CM filters, Muetze [48] used magnetic rings 
with relatively simple structure and convenient installation to suppress CM current. This 
method can suppress the circulating bearing currents and the rotor ground bearing cur- 
rents induced by the CM current but has a limited effect on the EDM bearing currents. 
In [15, 16], shielded motor cables are used to suppress rotor ground bearing currents. 

The second group comprises countermeasures within the motor. The most typical 
method is using insulated/ceramic/hybrid bearings to improve the impedance between 
the rotor and the stator frame, thus the bearing current coupling path of circulating bear- 
ing currents and rotor ground bearing currents are blocked [40]. However, considering 
the expense of ceramic bearings and the breakdown possibility of insulating bearings 
under high voltage and high frequency, large-sized motors are often equipped with 
shaft grounding brushes to eliminate the bearing currents [49]. In addition, a conductive 
ring mounted on the shaft which requires little maintenance compared to conventional 
grounding brushes is now widely spread in the industry [50]. It has been experimentally 
proven that bearing currents in the megahertz range such as EDM-bearing currents and 
circulating bearing currents can be effectively mitigated. In [51], the usage of conduc- 
tive grease is introduced to avoid bearing currents. However, the low-impedance path 
formed by the lubrication between the rolling elements and bearing race is not stable due 
to elastohydrodynamic effects. Besides, Busse [52] proposed an electrostatic shielded 
induction motor, which has an internally mounted Faraday shield to diminish the cou- 
pling between the stator and rotor. Experimental results show the modified motor solves 
the electrostatically induced bearing voltage without degrading its performance. Simi- 
larly, based on the electrostatic shield, grounded electrodes or slot wedges are introduced 
in stator slots to eliminate PWM-induced EDM bearing currents [53, 54]. 

Finally, different countermeasures and their effectiveness on different types of bear- 
ing currents are summarized in Table 1 (“represents effective, “x”represents not 
effective). 


Table 1. The effectiveness of different bearing currents mitigation techniques. 


Countermeasures Bearing current type 


EDM bearing Circulating bearing Rotor ground 
currents currents bearing currents 


Modified inverter f a/ A 
modulation strategies/CM 
voltage filters 


(continued) 
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Table 1. (continued) 


Countermeasures Bearing current type 
EDM bearing Circulating bearing Rotor ground 
currents currents bearing currents 
CM chokes/dv/dt filters x wh J 
Electrostatic shield/static afi x x 
charge 
dissipation/conductive 
grease 
Shielded cables/insulted x x xf 
couplings 
One x A x 
insulated/ceramic/hybrid 
bearing 
Two a g J 
insulated/ceramic/hybrid 
bearings 
One shaft grounding brush / x af 
Two shaft grounding J hr af 
brushes 
One J i J 


insulated/ceramic/hybrid 
bearing and one brush 


6 Conclusion 


In this paper, an overview of the current domestic and foreign researches on bearing 
currents in variable speed drive systems is presented. Specifically, the generation mech- 
anisms and equivalent circuit models for different types of bearing currents are sum- 
marized, and the theoretical calculation methods of capacitances related to the bearing 
current prediction are introduced. Finally, the commonly used mitigation techniques for 
bearing currents are discussed. Although considerable efforts have been made in this 


field, researches on relevant directions enumerated below are still urgently needed: 


(1) A universal equivalent circuit model for bearing current predictions. In this paper, 
the most typical four modes of bearing currents are discussed individually. While 
in practical applications, the topology of the overall inverter-motor system can 
vary greatly. Taking the metro vehicle traction system as an example, the negative 
dc-bus rail of the inverter and the motor frame are not directly grounded, which 
makes the modeling and prediction of the bearing currents particularly complicated. 
More work is needed to quantify the bearing current phenomena in these kinds of 
sophisticated systems. 
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(2) More low-cost and flexible mitigation techniques for bearing currents. In Part 5, 


various bearing current suppression methods have been summarized, while it is 
still a great challenge to turn these methods into industrial applications consid- 
ering the expense and reliability. Recent researches on new lubrication materials, 
which not only provide excellent lubrication characteristic but also suppress arc 
discharges with their high conductivity, are regarded as a solution to eliminating 
bearing currents. 
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Abstract. Polymer film capacitors are an efficient energy storage and conversion 
device, which has a wide range of applications in the field of electrical engineering. 
In this paper, a commercial BOPP film is selected as the dielectric film, and three 
metal electrode materials of Al, Cu, and Pt are grown on the surface of the BOPP 
film by vacuum evaporation or magnetron sputtering to explore the influence of 
metal electrode materials on the dielectric energy storage characteristics of BOPP 
films. This study found that compared with Al electrode films, BOPP films grown 
on Pt and Cu electrodes have lower leakage current density and higher breakdown 
field strength under high temperature conditions. When the applied electric field is 
630kV/mm, the maximum discharge energy density of Pt/BOPP/Pt film capacitors 
is 5.48 J/cm? at 125 °C, which is 1.34 times that of Al/BOPP/AI film capacitors, 
and the discharge efficiency remains above 90%. The results in this paper show that 
electrode materials with different metal work functions have significant effects on 
the dielectric properties of BOPP films. By choosing metal electrode materials with 
higher work functions can reduce the high-temperature leakage current density of 
BOPP films and increase their breakdown strength and energy storage density. 


Keywords: BOPP film - Dielectric properties - Metal electrode - Energy storage 
performance 


1 Introduction 


Film capacitors have a wide range of applications in the fields of electrical engineering 
and power electronics, such as filtering, voltage equalization, and energy storage [1]. 
The ability to release stored energy and generate large currents in a very short period of 
time has important applications in the pulsed power such as electromagnetic ejection. Of 
course, with the increasing power density of power equipment and electronic devices, 
the performance of film capacitors has put forward higher requirements, especially under 
extreme conditions such as high temperature and high electric field [2, 3]. For exam- 
ple, the commercial bi-oriented polypropylene (BOPP) film capacitors can withstand 
a maximum operating temperature of 105 °C. In order to meet the operating require- 
ments of DC support capacitors for new energy vehicle inverters, the additional devices 
must be added for cooling, which undoubtedly increases the system cost and physical 
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space. In addition, under the conditions of high temperature and strong electric field, 
the conductive loss of dielectric film increases sharply, leading to a serious deterioration 
of the discharge energy density (Ue) and discharge efficiency (7) of the capacitors [4, 
5], and the lost energy further transform into Joule heat, which accelerates the aging of 
dielectric film and even insulation failure, brings great hidden danger to the safe and 
stable operation of capacitors. 

Under high temperature conditions, Schottky charge injection is one of the main 
factors leading to the increased conductivity loss in thin film capacitors. Under the 
coupling conditions of high temperature and high electric field, the charges gain enough 
energy to overcome the interfacial barrier at the metal electrode and the dielectric film, 
and are injected into the interior of the dielectric [6-9]. The Schottky emission can be 
expressed by equation (1): 


Js = AT? exp (° + Jare VE) /koT (1) 


where EF is the applied electric field, A is the Richardson constant, T is the temperature. 
is the height of the potential barrier at the electrode/dielectric interface, € is the dielectric 
constant, and ko is the Boltzmann constant. The injection current density of the dielec- 
tric depends mainly on the height of the interfacial barrier, temperature, and the applied 
electric field. Polymers filled with wide bandgap nanofillers, such as silica (SiO2), alu- 
mina (Al203) and boron nitride (BN) nanosheets, introduce deep traps, limiting charges 
transport [10-12]. By growing an insulating or semiconducting layer with wide bandgap 
to polymer surfaces to form sandwich structures has proven to be an effective way to 
increase the height of the interfacial barrier at the electrode/dielectric. Li et al. deposited 
high bandgap width SiO2 on the surface of BOPP by CVD and PECVD, which hindered 
the charge injection from the electrodes and greatly reduced the space charge density, 
making BOPP film still possess high discharge energy density and efficiency at high 
temperatures [13]. Chi et al. introduced BN on both sides of polycarbonate (PC) and 
constructed a sandwich structure of BN/PC/BN, which effectively reduced the conduc- 
tion loss of polymer films at high temperatures and high fields [14]. Recently, Dong 
et al. compose a polymer dielectric with a laminar structure of wide bandgap inorganic 
Al203 and polyimide layers by solution reaction method, which still had a high discharge 
energy density of 1.59 J/cm? at 200 °C and discharge efficiency greater than 90% [15]. 
In addition, both phase field simulations and experimental results have shown that the 
introduction of a wide bandgap layer on the surface or inside of polymer can effectively 
inhibit charge injection and hinder the growth and diffusion of electric trees [16]. Under 
high temperature and high field, to enhance the height of the interfacial barrier between 
the polymer and the electrode is an effective way to achieve large energy storage density 
and high discharge efficiency. 

Film capacitors are mainly composed of polymer films and metal electrodes. The 
commonly used metal electrode materials are platinum (Pt), gold (Au), aluminum (Al), 
copper (Cu), zinc (Zn), etc. in laboratory, and their work functions are 5.61 eV, 4.32 
eV, 4.28 eV, 4.33 eV, and 4.33 eV, respectively. Metal electrodes with different work 
functions will have effects on the height of the interfacial potential barrier between the 
polymer and the electrode, but few related studies have been reported. Based on the above 
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mentioned, BOPP is selected as the dielectric material and Cu, Al and Pt are selected as 
the metal electrodes, respectively, in this study. Metal electrodes are grown on the upper 
and lower surfaces of BOPP films by vacuum vapor deposition or magnetron sputtering 
methods, including Al/BOPP/AI, Cu/BOPP/Cu, and Pt/BOPP/Pt capacitive films. The 
microstructure of BOPP films with different metal electrodes was investigated, and the 
mechanisms of the influence of different electrode materials on the dielectric, breakdown 
and energy storage properties of BOPP films at high temperatures have been clarified. 


2 Experimental materials and methods 


2.1 Experimental materials 


BOPP (thickness 5 um) film was purchased from Poly Technologies, USA; copper metal 
target was purchased from Hefei Kejing Material Technology Co. 


2.2 Preparation of BOPP film metallized electrodes 


Vacuum evaporation coating technology was used to grow Al on the upper and lower 
surface of BOPP film. Metallic Cu and Pt were grown on the top and bottom surfaces of 
the commercial BOPP films by magnetron sputtering technology. The oil and impurities 
on the surface of BOPP film were removed with anhydrous ethanol and deionized water. 
BOPP film was cut into small squares of 1 cm x 1 cm size and fixed on the mold and 
put into the high vacuum evaporation coating machine (ZHD-40 Beijing Techno Tech) 
for evaporation coating. High vacuum magnetron sputtering coater (JCP-350 Beijing 
Teknologue Technology Co., Ltd.) was employed to fabricate Cu and Pt electrodes, and 
the sputtering power was selected as 10 W and the sputtering time was 20 min (Fig. 1). 


(a) Vacuum vapor deposition of Al electrode (b) Magnetron sputtering of Pt and Cu electrodes 


Fig. 1. Schematic diagram of the preparation of films with different electrodes 


2.3 Microstructure characterization 


Scanning electron microscopy (SEM, SU8020, Japan) was used to characterize the 
microstructure of BOPP films and analyze their cross-sectional morphology. Atomic 
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force microscopy (AFM, CSPM, China) was used to analyze the surface roughness of 
different electrodes. X-ray diffractometer (XRD, Empyrean, Holland) was used to ana- 
lyze the phase structure of BOPP films. An infrared spectrometer (FTIR, Equinox55, 
Germany) was used to characterize the chemical bonding and molecular structure of the 
film samples. Differential scanning thermal analysis (DSC, Mettler Toledo, Switzerland) 
was used to analyze the thermal properties of BOPP films. 


2.4 Electrical performance characterization 


The films were tested in the frequency range from | Hz to 1 MHz using a broadband 
dielectric spectrum analyzer (Alpha-A, Novocontrol, Germany). The hysteresis curve 
(D-E curve) and leakage current density (I-V curve) of the films were tested using a fer- 
roelectric synthesizer (Precision Premier I, Radiant, USA). DC breakdown test system 
was used to test the breakdown strength of the films, and the accumulated breakdown 
probability was calculated by combining the Weibull distribution. 


3 Results and Discussion 


3.1 Microstructure 


The results of microstructural characterization of BOPP films are shown in Figure 2. 
From Figure 2(a), it can be seen that three sets of strong infrared absorption peaks appear 
in the wavelength range of 3500-500 cm7! in the infrared spectrum, among which the 
absorption peaks in the range of 1100-800 cm™! belong to the absorption band of fully 
homogeneous polypropylene crystals, the absorption peaks in the range of 1500-1300 
cm! belong to the C-H bending vibration mode of BOPP, and the range of 3000-2800 
cm! belongs to the C-H stretching vibration mode of BOPP [17]. Figure 2(b) shows the 
XRD diffraction pattern, from which it can be seen that the diffraction peaks at 14.1°, 
16.9°, 18.5°, and 25.5° in the diffraction angle range from 10° to 50° correspond to 
the (110), (040), (130), and (060) crystal planes, respectively, and the diffraction peak 
characteristics indicate that the BOPP films are of a-crystalline type [18]. As can be 
seen in Figure 2(c), the melting temperature of BOPP film is about 166.5 °C, and no 
obvious melting behavior occurs at 80-160 °C. On the one hand, this indicates that 
BOPP is highly crystalline, and on the other hand, it predicts that BOPP film has better 
structural thermal stability performance. Figure 2(d) shows the SEM structure of the 
cross-sectional BOPP film. As can be seen from the figure, the BOPP film section is flat 
and free of obvious defects, and it can also be observed that the thickness of the BOPP 
film is about 5 pm. 

Figure 3 shows the AFM images of BOPP films and BOPP films with different 
electrodes. It can be seen that the surface roughness of BOPP film is 7.32 nm, which 
of Al, Cu and Pt electrodes are 15 nm, 9.8 nm and 4.8 nm, respectively. The surface 
of Pt electrode is the lowest, where Al electrode is much larger than that of Cu and 
Pt electrodes, causing by the two different processes of vacuum vapor deposition and 
magnetron sputtering. 
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Fig. 3. AFM images of BOPP film and films with different electrodes 


3.2 Dielectric properties 


Figure 4 shows the dielectric frequency spectra of the BOPP films grown with dif- 
ferent electrodes. The results show that the relative permittivity of Pt/BOPP/Pt and 
Cu/BOPP/Cu capacitive films is slightly higher than that of Al/BOPP/AI films, while 
the dielectric loss of the former one is higher than the latter. The reason may be due 
to the preparation method of metallic Pt and Cu electrodes by magnetron sputtering. 
During sputtering, argon ions continuously bombard the metal target, and the resulting 
high-energy metal particles are continuously deposited on the film surface, which may 
lead to the phenomenon of broken bonds of molecular chains or generation of reactive 


Effect of electrode materials on dielectric properties 689 


radicals on the BOPP film surface. Under the alternating electric field, the reactive radi- 
cals enhance the polarization behavior of the film, and also increase the dielectric loss of 
the film. In contrast, Al electrode is prepared by using vacuum vapor deposition method. 
During the vapor deposition process, the Al melts and becomes vapor deposited on the 
surface of BOPP film, which has low energy for metal vaporization molecules and less 
influence on the surface structure of BOPP film. Therefore, the film of vaporized Al 
electrode exhibits low dielectric loss. 
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Fig. 4. Dielectric spectrum of BOPP films with different electrodes 


3.3 Breakdown Characterization 


DC voltages with different magnitudes were applied to the samples, and the effective 
values of the breakdown strength of the film samples were recorded. Based on the 
following equation: 


P(E) = 1 — exp(—(E/a)?) (2) 


The Weibull distribution and the Weibull breakdown field strength are obtained, 
as shown in Figure 5. Pt/BOPP/Pt film has the highest breakdown field strength of 
597.9 kV/mm, followed by Cu/BOPP/Cu film with 557.2 kV/mm, and the lowest is 
Al/BOPP/AI film with 534.3 kV/mm. Considering that the work function of the three 
electrodes gradually decreases (Pt > Cu > Al), the corresponding interfacial barrier 
height also decreases. Because of the higher work function of Pt, the electrons need 
to overcome the higher barrier to escape from the metal and inject into the dielectric, 
according to Equation (1), it is known that increasing the interfacial barrier height can 
effectively reduce the charges injection and improve the breakdown strength. 


3.4 Leakage current density Characterization 


Figure 6 shows that the leakage current density of the films increases with the increase of 
the applied electric field. At high temperature, the leakage current density of Pt/BOPP/Pt 
films is lower than that of Cu/BOPP/Cu and Al/BOPP/AI films, and when the applied 
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Fig. 5. Breakdown characteristics of BOPP films with different electrodes 


electric field is 200 kV/mm, the leakage current densities of the three electrode films (Al, 
Cu, and Pt) are 1.55 x 1077, 1.02 x 1077, and 6.71 x 1078 A/cm?, respectively, which 
is consistent with the variation law of Weibull breakdown field strength in Fig. 5(b). This 
phenomenon is mainly due to the metal electrode work function decision. In addition, 
the surface roughness of the metal electrode on the film surface also has some influence 
on the leakage density, as in Figure 3 the AFM results show that the Pt electrode has the 
lowest surface roughness and lower leakage current density. 
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Fig. 6. Leakage current density of BOPP films with different electrodes 


3.5 Polarization and energy storage performance 


Figure 7(a) shows the Schematic D-E hysteresis loop of the BOPP films with three 
different metal electrodes at 500 kV/mm and 125 °C. From Fig. 7(a), it can be seen 
that the polarization of film with different electrodes (Al, Cu, and Pt) increases with 
the increase of electric field, and the maximum polarization reaches 1.17, 1.37 and 
1.41 pC/cm?, respectively. Figure 7(b) shows the discharge energy density and energy 
storage efficiency of the BOPP films with different electrodes versus the applied electric 
field at 125 °C. The discharging efficiency of Pt/BOPP/Pt films is higher than that of 
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Cu/BOPP/Cu and Al/BOPP/AI films after the electric field exceeds 220 kV/mm, and 
the discharging efficiency of the films with Cu and Al electrode decreases more rapidly 
with the increase of field and gradually becomes less than 90%. Due to the influence of 
dielectric constant, breakdown field strength and dielectric loss, the Pt/BOPP/Pt film has 
the maximum discharge energy density under the same electric field, and the maximum 
discharge energy density can reach 5.48 J/cm? under 630 kV/mm, which is 1.34 times 
of that of Al/BOPP/AI, while the discharge efficiency remains above 90%. 
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Fig. 7. Polarization characteristics and energy storage performance of BOPP films with different 
electrodes 


4 Conclusion 


In this paper, metallic Al, Cu, and Pt electrode materials are successfully grown on the 
surface of BOPP film by vacuum vapor deposition and magnetron sputtering methods, 
respectively. Comparing the AFM surface structures of the films with different elec- 
trodes, it is found that the surface grains and roughness of the Al and Cu electrodes 
are relatively large, while the Pt electrode surface roughness is lower and the grains are 
fine. The electrical properties show that the relative permittivity of the BOPP films is 
about 2.25, and the permittivity of the Pt/BOPP/Pt films is slightly higher than that of 
the Cu/BOPP/Cu and Al/BOPP/AI films. Pt/BOPP/Pt film has the largest polarization 
at 125°C and 500 kV/mm, and the discharging efficiency is as high as 92%, and the 
maximum discharge energy density is 5.48 J/cm. 
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Abstract. Under the influence of the skin effect and the electromagnetic cou- 
pling between multiple conductors inside the cable, the distribution parameters of 
the cable shows frequency variation and high non-linearity, which makes it more 
difficult to describe the broadband impedance. The current impedance model of 
direct drive permanent magnet synchronous generator does not consider the effects 
of submarine cables impedance. In this paper, we (1) establish the impedance 
model of grid side converter, (2) extract distribution parameters of cable based 
on the COMSOL-MATLAB co-simulation platform, and (3) analyze the influ- 
ence of cable distribution parameters on wideband impedance characteristics of 
the wind power converter. The results show that the electromagnetic coupling 
between conductors mainly affects the cable distributed inductance and enhances 
the negative damping characteristics. After considering the distribution parame- 
ters, the submarine cable impedance presents multi-resonance characteristics. The 
longer the length is, the more obvious the multi-resonance characteristic will be, 
and the worse the broadband impedance characteristic of wind power converters 
connected to weak power grids will be. 


Keywords: Submarine cable - Distribution parameter - COMSOL-MATLAB - 
Wind power converter - Wideband impedance 


1 Introduction 


With the rapid development, more and more direct drive permanent magnet synchronous 
generators (D-PMSG) are connected to offshore wind farms. The submarine cable is an 
important part of the power collection network system, connecting wind turbines and the 
point of common coupling (PCC). When analyzing the stability of the system connected 
into the weak grid, the interaction between the submarine cable and the wind power 
system as a whole is usually observed as part of the weak power grid [1]. In modeling, 
the submarine cable is usually built as a resistance inductive element, and the distributed 
capacitance of the submarine cable is large, so the interactive coupling between the D- 
PMSG impedance and the grid impedance at a broadband frequency will stimulate more 
complex resonance characteristics. 
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The lumped parameter model is usually used to model transmission lines. Ref. [2] 
used PI-equivalent circuit to analyze the subsynchronous oscillation of wind farm. But 
there is a problem of too many PI sections for long transmission lines, which reduces 
the simulation efficiency [3]. At present, the oscillation frequency of wind farms ranges 
from 10 to 2,000 Hz [4, 5]. Influenced by skin effect, cable parameters have obvious 
frequency variation effect, and multi-resonance characteristics, which should be reflected 
by distributed parameter model [6]. 

Generally, parameters of single-core coaxial cable can be calculated by using the 
fitting formula [7]. But multi-conductors inside three-core cable are coupled with each 
other, so the formula cannot be directly applied [8]. At present, the transmission char- 
acteristics of multi-core armored cables are mainly obtained by measurement, which 
is not good for practical engineering application [9]. Therefore, a parameter extraction 
method suitable for multi-core submarine cable is needed. 

In this paper, the impedance model of D-PMSG GSC with submarine cable is estab- 
lished. Based on COMSOL-MATLAB co-simulation platform, the frequency variation 
distribution parameters of three-core submarine cable considering the coupling between 
conductors were extracted. The impedance transformation of lumped parameter model 
and distributed parameter model is compared by impedance model, and the effect of 
interphase coupling and cable length on impedance is analyzed. 


2 Modeling of Converter Considering Submarine Cable 


D-PMSG grid-connected system is shown in Fig. 1. It consists of wind turbine, PMSG, 
machine side converter (MSC), grid side converter (GSC), transformer and submarine 
cable. Among them, MSC controls the maximum power tracking, and GSC controls the 
voltage stability of the DC bus while adjusting the active/reactive power. 


Fig. 1. Diagram of grid connected D-PMSG 


Combined with harmonic linearization and symmetrical components, the balanced 
three-phase system can be decomposed into independent positive and negative order 
subsystems. When both subsystems are stable, the original three-phase system is also 
stable. The grid-connected system studied in this paper is three-phase symmetric, and 
the coupling effect between positive and negative sequence can be ignored. 


2.1 Impedance Model of D-PMSG 


Due to the existence of dc bus voltage regulator, the controlled current source can be 
used to equivalent to the change of wind turbine output power during stable operation, 


Impedance Modeling of Wind Turbine Converter Considering 695 


and the DC voltage is simplified to a constant value while ignoring the influence of the 
generator side. This paper is mainly based on the direct drive fan network side converter 
impedance modeling. 

The main circuit topology of the GSC is shown in Fig. 2. Where Uac is the DC 
voltage; Cac and Rac are the voltage regulator capacitor and resistance; Ky is the voltage 
control transfer function; Gj is the current control transfer function; H; is the equivalent 
transfer function of current sampling delay, zero order hold and low pass filter. Ref. [10] 
resumes the positive sequence and negative sequence impedance of the GSC considering 
the influence of the PLL based on the harmonic linearization: 
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Fig. 2. Diagram of grid connected D-PMSG converter 
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where, I1 = Ia _ref/2 + jlq_ret/2. Ia_ref and Iq ref are the reference current of the power 
loop. 


2.2 Impedance Model of Transformer 


Transformer has the function of impedance transformation. According to the concept of 
input impedance, the fan impedance is converted to the t-type equivalent circuit of the 
primary side rear transformer, as shown in Fig. 3. 

The input impedance frequency domain expression is: 


[kG (R2 + sL2 + Zesc)|(Rm + sLm) 


(3) 
kå (Ra + sL + Zasc) + Rm + sLm 


Zra = Ri + sly + 


Where kr is the transformer ratio and Zgsc is the impedance of GSC. 
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Fig. 3. T type equivalent circuit of transformer 


2.3 Impedance Model of Submarine Cable 


At high frequency, distribution parameters of cable should be considered. The cable can 
be thought of as consisting of many infinitesimal dx (Fig. 4). 
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Fig. 4. Distribute model of submarine cable 


k and m represent the head and end of the line respectively. ro, lo, co, and go represent 
resistance, inductance, capacitance, and conductance per unit length, respectively. For 
a well-insulated cable, the leakage current is very small, so the shunt conductance can 
be ignored. The relationship between the voltage and current at the head and end of the 
cable in the frequency domain is: 
| Ùk | _ | cosh(y (w)L) Ze (œ) sinh (y (w)L) | d (4) 
ht ey ohyo) YL din 
Where: y(@) = /Zo(@)Yo(@) is propagation coefficient; Z.(w) = /Zo(@)/Yo(@) is 
characteristic impedance; L is the length of cable; Zo(@) = ro(w) + j@lo(@) is unit 
impedance; Yo(w) = jwCo(w) is unit admittance. In combination with Eqs. (3) and (4), 
the impedance model of D-PMSG with submarine cable distribution parameters is: 


_ Z[@re + Ze) + (Zra — Ze)A?] 


BS (Zro + Ze) — (Zro — Ze)A? (3) 


3 Distribution Parameters of Submarine Cable 


When electrical simulation software such as PSCAD and Simulink uses the line constant 
program (LCP) to compile the three-core cable, usually ignore the phase coupling. There 
is electromagnetic coupling between the multiple conductors in the multi-core submarine 
cable, and the skin effect is more complicated at high frequency. Ignore the coupling 
between phases to make the calculation error larger. 
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The finite element method is suitable for the calculation of cable parameters of 
various cross-sectional structures, and is not affected by the proximity effect, and is 
more accurate for the calculation of the incremental impedance at high frequencies. This 
chapter extracts the cable distribution parameters based on the finite element method, 
and studies the influence of the phase coupling on the cable impedance. 


3.1 Finite Element Model of Submarine Cable 


Common three-core armored submarine cables include a total of 7 metal structures: 
three-phase cable core, three-phase metal shielding layer and armor. Based on the elec- 
tromagnetic field theory and the 2D finite element model, the parameters can be easily 
obtained. The result of grid dissection is shown in Fig. 5. 


Fig. 5. Subdivision diagram of three-core armoured marine cable 


3.2 Multi-conductor Distributed Impedance Extraction 


Apply Ix to conductor i, according to the definition of resistance and inductance: 


_ Pa _ ff dsEJ? _ 1 Sf aside? 
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The numerator integration domain in Eq. (6) is all conductors, and the denominator 
integration domain is the conductor excited. Equation (7) considers the coupling between 
multiple conductors, and the molecular integration domain is for all conductors. 

Applying a forward current J, to conductor i and a reverse current - Iy to conductor j, 
the mutual impedance is obtained by the loss in the domain P4 and the magnetic energy 
W n, Le. 


1 
Rij = Ri = 5 (Ri +Rj— 4Pa IR) 
i (7) 
Ly = Li = 5 (Li + Li — 8Wn/I2) 
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The COMSOL LiveLink for MATLAB function is used to extract the multi-conductor 
system distribution parameters automatically. 


3.3 Submarine Cable Sequence Impedance Parameters 


The distribution parameters R, L, and C extracted based on the calculation method 
proposed above are all 7th order matrices. According to the series voltage drop on the 
cable core, metal shield and armor layer can be obtained from the full-rank matrix as 


Vo Ic 
vs |= Zee) ZCS,CM(3x4) I (8) 
Vu Z¢s,cm(3x4) ZSS,MM4x4 || 7 


where Vc = [Va Vp Vcl!, Vs = [Vsa Vss Vsc]!, Ic = Ua Is Ic", Is = Usa Isp 
Isc]. 


Zphase = Zcc — Zcs,smZ5. um Zsc,sM (9) 


where Zphase is a3 x 3 order core triple equivalent impedance matrix. 


4 Case Study 


In this paper, 5 MW direct-drive wind turbine is used as the research object to compare 
and analyze the effect of transformer as well as submarine cable on the impedance of 
direct-drive wind turbine. The influence of interphase coupling of submarine cables 
on the distribution parameters is considered, and the influence of interphase coupling 
and different lengths of submarine cables on the impedance characteristics of the grid- 
connected system of direct-drive wind turbines is compared and considered or ignored. 


4.1 Analysis of Submarine Cable Distribution Parameters 


Figure 6 compares different models, including PI-equivalent circuit, fixed distributed 
parameter, and frequency-dependent parameter model. The cable length is 10 km. 


Amp/dB 


Fig. 6. Comparison of impedance of D-PMSG with different cable models 


The impedance models in Fig. 6 are converted to the 35 kV side. A Fig. 6 section of 
the total parameter PI- equivalent circuit is used to equate 10 km cable, which has only 
one resonant peak at about 1 kHz, and the system impedance changes from inductive 
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to capacitive at about 25 Hz. The impedance model with fixed distributed parameter 
model in Fig. 6 has four resonance peaks of equal amplitude in 0 ~ 4 kHz, and the first 
resonance peak appears at a similar position to the PI- equivalent model. 

When considering the distribution parameter frequency change characteristics, the 
parameter change effect is not obvious at lower frequency, and the distribution parame- 
ter model curve in Fig. 6 has higher overlap in the low frequency band, and the system 
impedance changes from inductive to capacitive at about 15 Hz. At higher frequen- 
cies, the parameter change is obvious, and the multi-resonance peak considering the 
frequency-dependent effect appears later and the amplitude decays gradually. 


4.2 Analysis of Phase Coupling 


The parameters change significantly when frequency is high. To observe the skin effect, 
sweep the cable finite element model from 0 to 4 kHz. In order to compare the inter- 
phase coupling effect on the impedance of the grid-connected system, the impedance 
model with and without coupling is compared, and the distributed parameter model with 
frequency variation is used for the cable, and the lengths are set to 1 km and 5 km. the 
impedance model results are shown in Fig. 7. 


Amp/dB 


Fig. 7. D-PMSG impedance of submarine cable interphase coupling 


In Fig. 7, the negative damping characteristics of the cable are enhanced after the 
interphase coupling is considered, and the overall impedance of the high-frequency sys- 
tem is reduced. The low frequency band phase margin safety zone is reduced. The second 
resonance peak of the 10 km cable at greater than 1 kHz is due to the resonant charac- 
teristics of the cable itself. As the length increases, the influence of the multi-resonance 
characteristics gradually becomes obvious, and the broadband stability decreases when 
connected to a weak power grid. 


4.3 Analysis of Cable Length 


It is known from the previous section that the multi-resonance characteristics of long 
distance submarine cable under wide frequency cannot be neglected. In order to study 
the influence of the length of the submarine cable on the impedance characteristics of 
the grid-connected subsystem, this chapter adopts the distributed parameter submarine 
cable model considering the frequency variation effect and inter-phase coupling, and 
sets the lengths as 1, 5, 10 and 20 km, respectively. 

The effect of the submarine cable on the system impedance in Fig. 8 becomes more 
and more obvious with increasing length. The safety domain of the low frequency band 
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Fig. 8. Impedance of D-PMSG with different length of marine cable 


of the long-distance cable decreases, and for a grid-connected wind turbine system 
connected to a cable 20 km long, the system impedance changes abruptly to capacitive 
at frequencies above about 10 Hz. The resonance peaks in the high frequency band 
of long-distance submarine cable appear earlier and more. 1 km submarine cable has 
no obvious resonance peaks in the range of 4 kHz, 5 km submarine cable has the first 
resonance peak at about 3 kHz. 20 km submarine cable has the first resonance peak at 
about 700 Hz, and there are 5 obvious resonance peaks in 0—4 kHz, and the amplitude 
gradually decays. 


5 Conclusion 


In this paper, a direct-drive wind turbine grid-side converter impedance model with 
cable distribution parameters is established, and a three-core submarine cable finite 
element model is established. driven wind turbine grid-connected system by comparing 
the impedance model simulation results. The main conclusions are as follows. 


(1) The longer the length of the cable, negative damping and multi-resonant character- 
istics at high frequency more significantly. Taking into account the influence of the 
submarine cable, the system impedance increases capacitively at high frequencies, 
and the amplitude and phase angle are significantly reduced. 

(2) Three-core submarine cable inter-phase coupling effect mainly affects the subma- 
rine cable distribution inductance. In this paper, the joint COMSOL-MATLAB 
simulation model can be used to quickly extract the distribution parameter matrix 
of three-core armored submarine cable considering conductor coupling and skin 
effect. 

(3) A grid-connected impedance model of direct-drive wind turbines with subma- 
rine cable distribution parameters is established, which can be used to analyze 
the broadband stability between wind farms and the grid. 
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